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Abstract 



We consider QCD radiative corrections to the production of W and Z bosons in hadron 
collisions. We present a fully exclusive calculation up to next-to-next-to-leading order 
(NNLO) in QCD perturbation theory To perform this NNLO computation, we use 
a recently proposed version of the subtraction formalism. The calculation includes 
the j—Z interference, finite-width effects, the leptonic decay of the vector bosons and 
the corresponding spin correlations. Our calculation is implemented in a parton level 
Monte Carlo program. The program allows the user to apply arbitrary kinematical 
cuts on the final-state leptons and the associated jet activity, and to compute the 
corresponding distributions in the form of bin histograms. We show selected numerical 
results at the Tevatron and the LHC. 
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The production of W and Z bosons in hadron collisions through the Drell-Yan (DY) mechanism 
[I] is extremely important for physics studies at hadron colliders. These processes have large 
production rates and offer clean experimental signatures, given the presence of at least one high- 
er lepton in the final state. Studies of the production of W bosons at the Tevatron lead to precise 
determinations of the W mass and width [2]. The DY process is also expected to provide standard 
candles for detector calibration during the first stage of the LHC running. 

Because of the above reasons, it is essential to have accurate theoretical predictions for the 
vector-boson production cross sections and the associated distributions. Theoretical predictions 
with high precisions demand detailed computations of radiative corrections. The QCD corrections 
to the total cross section [3] and to the rapidity distribution [I] of the vector boson are known up 
to the next-to- next-to-leading order (NNLO) in the strong coupling as- The fully exclusive NNLO 
calculation, including the leptonic decay of the vector boson, has been completed more recently 
[5]. Full electroweak corrections at 0(a) have been computed for both W [6] and Z production [TJ. 

In this Letter we present a new computation of the NNLO QCD corrections to vector boson 
production in hadron collisions. The calculation includes the j-Z interference, finite-width effects, 
the leptonic decay of the vector bosons and the corresponding spin correlations. Our calculation 
parallels the one recently completed for Higgs boson production [8j [9], and it is performed by 
using the same method. 

The evaluation of higher-order QCD corrections to hard-scattering processes is complicated by 
the presence of infrared (IR) singularities at intermediate stages of the calculation that prevents a 
straightforward implementation of numerical techniques. Despite this difficulty, general methods 
have been developed in the last two decades, which allow us to handle and cancel IR singularities 
[TUt [TT| [12] appearing in NLO QCD calculations. In the last few years, several research groups 
have been working on extensions of these methods to NNLO [TBI [HI [151 EE E], an d; recently, 
the NNLO calculation for e + e~ — > 3 jets was completed by two groups [HI [19]. Parallely, a new 
general method [20J, based on sector decomposition [21], has been proposed and applied to the 
NNLO calculations of e + e~ — *> 2 jets [22], Higgs [23] and vector [5] boson production in hadron 
collisions, and to some decay processes [21]. Our method [8] applies to the production of colourless 
high-mass systems in hadron collisions, and is based on an extension of the subtraction formalism 
[ITj[T2] to NNLO that we briefly recall below. 

We consider the inclusive hard-scattering reaction 

h 1 + h 2 ^V(q) + X, (1) 

where the collision of the two hadrons hi and hi produces the vector boson V (V = Z/j*, W + or 
W~), with four-momentum q and high invariant mass \fq^. At next-to- leading order (NLO), two 
kinds of corrections contribute: i) real corrections, where one parton recoils against V; ii) one- 
loop virtual corrections to the leading order (LO) subprocess. Both contributions are separately 
IR divergent, but the divergences cancel in the sum. At NNLO, three kinds of corrections must 
be considered: i) double real contributions, where two partons recoil against V; ii) real-virtual 
corrections, where one parton recoils against V at one-loop order; iii) two-loop virtual corrections 
to the LO subprocess. The three contributions are still separately divergent, and the calculation 
has to be organized so as to explicitly achieve the cancellation of the IR divergences. 

We first note that, at LO, the transverse momentum of V is exactly zero. As a consequence, 
as long as 7^ 0, the (N)NLO contributions are actually given by the (N)LO contributions to 
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V + jet(s). Thus, we can write the cross section as 
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This means that, when qx 7^ 0, the IR divergences in our NNLO calculation are those in da^Q S : 
they can be treated by using available NLO methods to handle and cancel IR singularities (e.g., 
the general NLO methods in Refs. [TOj [TTJ [12]). The only remaining singularities of NNLO type 
are associated to the limit — > 0. Following Ref. [8J we treat them by an additional subtraction. 
Our key point is that the singular behaviour of da^^Q when qx — > is well known: it comes 
out in the resummation program [25J of logarithmically-enhanced contributions to transverse- 
momentum distributions. Therefore, the additional subtraction can be worked out by using a 
counterterm, da?£ LO , whose general structure [8] depends only on the flavour of the initial- 
state partons involved in the LO partonic subprocess (qq annihilation in the case of vector-boson 
production, gg fusion in the case of Higgs boson production). 



Our extension of Eq. (J2J to include the contribution at qx = is [8] : 
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we have subtracted the (N)LO counterterm da? 7 



Comparing with the right-hand side of Eq. wc nave ouuuc»,rcu ^ v ±x ;uu wuuociucnu ^(^lo 
and added a contribution at qx = 0, which is needed to obtain the correct total cross section. 
The coefficient 7~C£ N } NLO does not depend on q T and is obtained by the (N)NLO truncation of the 
hard-scattering perturbative function 
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According to Eq. ([3]), the NLO calculation of da v requires the knowledge of H v( -^ and the LO 
calculation of da v+iets . The general (process-independent) form of the coefficient Ti^ is known: 
the precise relation between US 1 ' and the IR finite part of the one-loop correction to a generic 
LO subprocess is explicitly derived in Ref. [2D]. At NNLO, the coefficient H v ^ is also needed, 
together with the NLO calculation of da v+ i cts . The calculation of the general structure of the 
coefficients H.^ is in progress. Meanwhile, by using the available analytical results at C?(a§) f° r 
the total cross section j3] and the transverse-momentum spectrum [27| of the vector boson, we have 



explicitly computed the coefficient 7i v ^ of the DY process. Since the NLO corrections daY +1 ' >[> 
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to qq — > V+jet(s) are also available [28], using Eq. 
NNLO calculation of vector-boson production. 



we are able to complete our fully-exclusive 



We have encoded our NNLO computation in a parton level Monte Carlo program, in which 
we can implement arbitrary IR safe cuts on the final-state leptons and the associated jet activity. 

In the following we present an illustrative selection of numerical results for Z and W production 
at the Tevatron and the LHC. We consider u, d, s, c, b quarks in the initial state. In the case of W ± 
production, we use the (unitarity constrained) CKM matrix elements V u d = 0.97419, V us = 0.2257, 
V ub = 0.00359, Vd = 0.2256, V cs = 0.97334, V cb = 0.0415 from the PDG 2008 [29]. In the case of Z 
production, additional Feynman diagrams with fermionic triangles should be taken into account. 
Their contribution cancels out for each isospin multiplet when massless quarks are considered. The 
effect of a finite top-quark mass in the third generation has been considered and found extremely 
small [3D], so it is neglected in our calculation. As for the electroweak couplings, we use the so 
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called scheme, where the input parameters are Gp , rn Z} m w . In particular we use the values 
G F = 1.16637 x 1(T 5 GeV" 2 , m z = 91.1876 GeV, T z = 2.4952 GeV, m w = 80.398 GeV and 
F w = 2.141 GeV. We use the MSTW2008 [3]] sets of parton distributions, with densities and a s 
evaluated at each corresponding order (i.e., we use (n + l)-loop as at N n LO, with n = 0, 1,2). 
The renormalization and factorization scales are fixed to the value /i# = fiF = rnvi where my is 
the mass of the vector boson. 

We start the presentation of our results by considering the inclusive production of e + e~ pairs 
from the decay of an on-shell Z boson at the LHC. In Fig. [1] (left panel) we show the rapidity 
distribution of the e + e~ pair at LO, NLO and NNLO, computed by using the MSTW2008 partons. 
The corresponding cross sections^ are a^o = 1.761 ± 0.001 nb, a NLO = 2.030 ± 0.001 nb and 
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Figure 1: Rapidity distribution of an on-shell Z boson at the LHC. Results obtained with the 
MSTW2008 set (left panel) are compared with those obtained with the MRST2004 set (right panel) . 

o~nnlo — 2.089 ± 0.003 nb. The total cross section is increased by about 3% in going from NLO 
to NNLO. In Fig. [T] (right panel) we also show the results obtained by using the MRST2002 LO 
[32] and MRST2004 [53] sets of parton distribution functions. The corresponding cross sections 
are a LO = 1.629 ± 0.001 nb, a NLO = 1.992 ± 0.001 nb and a NNLO = 1.954 ± 0.003 nb. In this case 
the total cross section is decreased by about 2% in going from NLO to NNLO. 

We next consider the production of e + e~ pairs from Z/'-f* bosons at the Tevatron. For each 
event, we classify the lepton transverse momenta according to their minimum and maximum 
values, prmin and prmax- The leptons are required to have a minimum p T of 20 GeV and pseudo- 
rapidity \r]\ < 2. Their invariant mass is required to be in the range 70 GeV < m e + e - < 110 GeV. 

^Throughout the paper, the errors on the values of the cross sections and the error bars in the plots refer to an 
estimate of the numerical errors in the Monte Carlo integration. 
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Figure 2: Distributions in prmin o,nd prmax for the Z signal at the Tevatron. 

The accepted cross sections are o^o = 103.37 ± 0.04 pb, (Jnlo = 140.43 ± 0.07 pb and <Jnnlo = 
143.86 ± 0.12 pb. In Fig. [2] we plot the distributions in prmm and prmax at LO, NLO and NNLO. 
We note that at LO the prmm and prmax distributions are kinematically bounded by pr < Q max/2, 
where Q max = 110 GeV is the maximum allowed invariant mass of the e + e~ pairs. The NNLO 
corrections have a visible impact on the shape of the prmm and prmax distribution and make the 
VTmin distribution softer, and the PTmax distribution harder. 

We finally consider the production of a charged lepton plus missing p? through the decay of a 
W boson (W = W + , W~) at the Tevatron. The charged lepton is selected to have pt > 20 GeV 
and \tj\ < 2 and the missing pt of the event should be larger than 25 GeV. We define the 
transverse mass of the event as tut = — cos0), where cf> is the angle between the 

the pt of the lepton and the missing px- The accepted cross sections are cr^o = 1.161 ± 0.001 nb, 
a NLO = 1.550 ± 0.001 nb and a NNLO = 1.586 ± 0.002 nb. 

In Fig. [3] we show the rriT distribution at LO, NLO and NNLO. We note that at LO the 
distribution has a kinematical boundary at my = 50 GeV. This is due to the fact that at LO 
the W is produced with zero transverse momentum: therefore, the requirement p^ lss > 25 GeV 
sets tut > 50 GeV. Around the region where = 50 GeV there are perturbative instabilities 
in going from LO to NLO and to NNLO. The origin of these perturbative instabilities is well 
known [31]: since the LO spectrum is kinematically bounded by my > 50 GeV, each higher- 
order perturbative contribution produces (integrable) logarithmic singularities in the vicinity of 
the boundary. We also note that, below the boundary, the NNLO corrections to the NLO result 
are large; for example, they are about +40% at mx ~ 30 GeV. This is not unexpected, since 
in this region of transverse masses, the 0(a$) result corresponds to the calculation at the first 
perturbative order and, therefore, our 0(a$) result is actually only a calculation at the NLO level 
of perturbative accuracy. 
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Figure 3: Transverse mass distribution for W production at the Tevatron. 



We have illustrated a calculation of the cross section for W and Z boson production up to 
NNLO in QCD perturbation theory. An analogous computation was presented in Ref. j5]. Our 
calculation is performed with a completely independent method. In the quantitative studies that 
we have carried out, the two computations give results in numerical agreement. Our calculation is 
directly implemented in a parton level event generator. This feature makes it particularly suitable 
for practical applications to the computation of distributions in the form of bin histograms. A 
public version of our program will be available in the near future. 

Acknowledgements. We thank Stefan Dittmaier for useful correspondence. 
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